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Abstract 

A high order linear (/-difference equation with polynomial coefficients 
having g-Hahn multiple orthogonal polynomials as eigenfunctions is given. 
The order of the equation is related to the number of orthogonality condi- 
tions that these polynomials satisfy. Some limiting situations when q 1 
are studied. Indeed, the difference equation for Hahn multiple orthogonal 
polynomials given in [7] is corrected and obtained as a limiting case. 

1 Introduction 

The relevance of the special functions and classical orthogonal polynomials as 
eigenfunctions of a second order differential equation is a well established fact 
[5]. For instance, many fundamental problems of quantum mechanics as the 
harmonic oscillator, the solution of the Schrodinger, Dirac and Klein-Gordon 
equations for a Coulomb potential, the motion of a particle in homogeneous 
electric or magnetic field lead to the study of the eigenfunctions for the gener- 
alized equation of hypergeometric type (see 0). Recall that the spherical and 
cylinder (Bessel) functions -perhaps the most popular special functions- as well 
as the classical orthogonal polynomials are particular solutions of this equation 
(see also equation (TT]) below). 

Recently, the appearance of new special functions, namely, the multiple or- 
thogonal polynomials also having the property of being eigenfunctions of a dif- 
ferential equation have attracted the interest of several researchers -see [5] , and 
,f; for a discrete case-. In an effort to consider more general situations, the 
notion of g-Hahn multiple orthogonal polynomials was introduced in [3]. Fur- 
thermore, a particular situation in which a third order g-difference equation 
having the aforementioned multiple orthogonal polynomials as eigenfunctions 
was considered. Here we complete the study initiated in [3] by obtaining a 
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high order g-difference equation having g-Hahn multiple orthogonal polynomi- 
als as eigenf unctions. This (/-difference equation gives a relation for the poly- 
nomial with a given degree evaluated at non- uniformed distributed points x(s), 
x(s + 1), . . . , x(s + r + 1); where x(s) denotes the g-exponential lattice and 
r, the number of orthogonality conditions In addition, other multiple 

orthogonal polynomial families like some studied in [H Q~2] as well as their cor- 
responding differential and difference equations (see [5JE]) can be obtained as 
a limiting case. 

The contents of this paper are as follows. In Section some well-known 
results on classical and multiple orthogonality are summarized. In Section [3] 
we give an explicit expression for the multiple orthogonal polynomials studied 
here as well as a high order g-difference equation that these polynomials satisfy. 
Indeed, the notion of lowering and raising operators as well as the auxiliary 
lemma [3~T1 play an essential role in the accomplishment of the g-difference equa- 
tion (theorem 13. ip . As corollary we get some previous results [7] . Lastly, 
Section |4] comprises some limiting cases known in the literature as multiple Ja- 
cobi and Hahn polynomials, respectively [H 03 [T^]. We also show how these 
families of polynomials can be obtained from the g-Hahn multiple orthogonal 
polynomials based on the Rodrigues-type formula. 



2 Preliminary notions 

One of the most remarkable features of the classical orthogonal polynomials on 
the real line, namely, Jacobi, Laguerre, and Hermite polynomials is that they 
are eigenfunctions of the second order differential operator [TJ [6] 

f d? d \ 

\ a2 ^dx^ + ai ^d^ +XI j y( - X ^ 0, (1) 

where deg 0,2(2;) < 2, degai(a;) = 1, A is a constant independent on x and X 
denotes the identity operator. The above equation is known as hypergeome- 
tric equation. An exhaustive study of their solutions in term of the polynomial 
coefficients 0.2(0;) and ai(a;) for each classical family can be found in [5]. 

The notion of orthogonality for these classical families is assumed to be with 
respect to a Borel measure \i on the real line M (with infinitely many points 
of increase) supported on a subset ft C K, where U — (a, b), with \a\ < 00, is 
the smallest interval on the real line that contains supp /x. This orthogonality 
concept can be briefly described as follows: If f n p(x)dn(x) converges for every 
polynomial p, then one can define an inner product 

(p,q) = [ p(t)q(t)d»(t), 

where p, q are polynomials. For such an inner product, which we will say to 
be standard, a sequence of polynomials (p n )n>o is said to be orthogonal with 
respect to the above inner product if 

(i) degp n (x) = n, 

(ii) (p n ,Pm) = dn.m lbn|| 2 , tti, u € N, and S m , n denotes the Kronecker delta 
function. 
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A difference analog of the equation (fTJ) on the lattice x(s), {x(s) i— » K + : 
a < s < b — 1}, is the hypergeometric-type difference equation [5] 

y(s) = y(x(s)), a(s) = a 2 {x(s)) - ±ax(x(s)) Ax(s- a), r(s) = oi(x(s)), 

(2) 

being V2/( s ) = Ay(s — 1). Here A and V are the forward and backward differ- 
ence operators, respectively. Observe that A/(s) = /(s + 1) — /(s). 

Analogously, the classical orthogonal polynomials of a discrete variable can 
be obtained as the corresponding eigenfunctions of ^j. Indeed, the polynomial 
family that verifies ^ can be orthogonalized by constructing a Sturm-Liouville 
problem with orthogonalizing weight w(s) as the solution of a Pearson-type 
difference equation (see [8j pp. 70-72]). Indeed, the following orthogonality 
property 

6-1 

J2Pn{x{s))P m (x(s))u(s) A x(s - |) = <5 n , m ||.P„|| 2 , (3) 

s—a 

yields under the additional (boundary) conditions a(s)uj(s)x k (s — j)| _ b = 0, 
k = 0, 1, . . . From now on we will denote any polynomial P n (x(s)) simply as 

Pn(s). 

In [5] was studied the g-Hahn orthogonal polynomials P"'°(s), n = 0, 1, . . . , 
with parameters a, (3 > —I, on the lattice 

x( s ) = ^l q ^l. (4) 

In particular, the orthogonalizing weight and the coefficients of the second order 
difference equation are as follows: 



oa ,fS,N, s) = T q (s + P+l)r q (N + a-s + l) 

u(*) = { q q-^ s r q ( s + i)r q (N + i- s ) 



00 = ,±k - , - — - « = o,i,...,JVeN, 

0, otherwise, 



(5) 

q < 1, or f (a) = /(s;^ 1 ) if ? > 1, being /(*;<?) = (1 - g )i-^ ^>o^^) and 



where N denotes the set of natural numbers, f(s) = q ' * 'f(s',q) if < 



<r(s) = -g N 2° x ( s ) 2 + g 2 [JV + a] g ar(s), where [x], — ' J '' 



j<? i _i j 
q? ~ q 2 

r( s ) = -< z i!± ^ [a + /3 + 2],i(«) + q a+fj+1 [P + l] q [N - l] q , (6) 

\ n = q * [n] q [n + a + (3 + l] q . 

For the classical orthogonal polynomials there is a characterization due to 
Sonin (and independently obtained by W. Hahn) in the sense that they are the 
families of orthogonal polynomials such that the sequence of their first deriva- 
tives constitutes again a sequence of orthogonal polynomials with a shifted set of 
parameters. For classical orthogonal polynomials of a discrete variable there is a 
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similar characterization. When one takes the forward difference of the classical 
orthogonal polynomials of a discrete variable, one can show that these polyno- 
mials are again orthogonal polynomials of the same family, but with a different 
set of parameters. As such, the Ax ^_ operator acts as a lowering operator 
on these families of polynomials. As an example for the g-Hahn polynomials 

~ r ,P^ N (s) = q -nn] q P:+l^ N -\s). (7) 



Ax(s 

Also for the classical orthogonal polynomials of a discrete variable we have 
raising operators which can also be found using the orthogonality property Q 
and summation by parts. Hence, 



1 V 



a,f},N (s) pa ,f3,N (s) 



a - 



q~ 



' Zll pa-l,P-l,N+l( \ 
v + n + g r n+l \ a J- 



(8) 



By combining the lowering and raising operators (JT])-© one can obtain the 
equation {3J for the g-Hahn orthogonal polynomials with the above polynomial 
coefficients ©. This fact is a general feature of the classical orthogonal poly- 
nomials of a discrete variable, i.e., the corresponding combination of lowering a 
raising operators leads to the second order difference equation (|2|). 

In section[3]we will generalize the above procedure for a non-standard orthog- 
onality. Indeed, a new situation appears, namely, a third order ^-difference equa- 
tion having g-orthogonal polynomials as eigenfunctions. This situation keeps 
certain similarity with the previous accomplishment of the hypergeometric equa- 
tion via raising and lowering operators; however two orthogonality conditions 
must be considered instead (see j3]). Hence, we need to deal with the concept of 
multiple orthogonal polynomial. This notion appears naturally in simultaneous 
Pade approximation -which is often known as Hermite-Pade approximation- in 
order to get simultaneous rational approximants to a vector function. As we 
will see below they can be interpreted as an extension of ordinary orthogonal 
polynomials. 

Let /ii, fj,2, ■ ■ ■ , fJ-r (f > 2) be Borel measures on R with infinitely many 
points of increase such that supp /jj C fl<, i = l,...,r, where 51^ = (aj, with 
a,; | < oo, is the smallest interval on the real line that contains supp/i^. The 
Cauchy transforms of the corresponding measures 

fi-i(z) = [ ^IhJL ; z^fti i = l,...,r, 



also known as Markov (or Stieltjes) functions, can be simultaneously approxi- 
mated by rational functions with prescribed order near infinity. For studying 
this problem [10) one needs to introduce a multi-index ft = (ni, ri2, • • ■ , n r ) 
of nonnegative integers, and finds a polynomial Pn(z) ^ of degree at most 
\n\ = n\ + • • • + n r , such that the expressions 

P*(z)pLi{z) - Qn,i(z) = + ■■■ = 0{z- n '- v ), i = 1, . . . , r, 
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hold. Observe that the polynomial Pft(z) is the common denominator of the 
simultaneous rational approximants Qn,i{ z )/ Pn( z )i to the Markov (Stieltjes) 
functions fli(z), i = 1, 2, . . . ,r. 

Indeed, Pft(z) is a type II multiple orthogonal polynomial of degree < \n\ 
defined by the following orthogonality conditions (see [4l ITOl [12] ) 

/ Pn{x)x k d^i{x) = 0, k = 0, l,...,m - 1, j=l,...,r. (9) 

These conditions Q give a linear system of |n| homogeneous equations for the 
|n| + 1 unknown coefheients of Pa(z). If the multi-index n is normal |10j the 
solution is a unique polynomial Pn(z) (up to a multiplicative factor) of degree 
exactly \n\. In this situation throughout the paper we consider always monic 
multiple orthogonal polynomials. 

If the measures in © are positive discrete measures on K, i.e., 

Ni 

Hi = / ]vi,k5xi, h , > 0, Xi : k £ K, iV l eNU{+oo}, i = l,...,r, 

where Xj x ^ 7^ :Ei 2l fc, A; = 0, . . . , A^, whenever i\ ^ the corresponding polyno- 
mial solution is then a discrete multiple orthogonal polynomial Pn{z) [3]. Here 
we have that supp/i; is the closure of {xi t k}k=o an( ^ that is the smallest 
closed interval on R that contains {xi t k}k=o- ^ the a bove system of measures 
forms an AT system then every multi-index is normal (see [10) for a detailed 
explanation on the concept). 

In [1] for several AT systems of measures was studied the corresponding 
discrete multiple orthogonal polynomials of type II on the linear lattice x(s) = s 
(those of Charlier, Kravchuk, Meixner of first and second kind, and Hahn). Even 
more, it was obtained rising operators and then the Rodrigues-type formula. 

Here we will deal with an AT system formed by different discrete measures 
supported on the same interval. This situation was analyzed in [3] for n = 
(n\, . . . , n r ) £ N r and the set of parameters N, a>o and a — (ct\, . . . , a r ). Indeed, 
the orthogonality conditions ^ were considered with respect to the following 
positive discrete measures on K 

N 

(ii = J2v« i ' ao ' N (k)Ax{k-h)S(k-s), i = l,...,r, (10) 

s=0 

where v^ i ' a °' N (k) is defined in ([5]) and cto, on > —1, cti — a>j £ {0, 1, . . . , N — 1} 
when i ^ j. 

Definition 2.1. A polynomial p"' a °> (s) that verifies the orthogonality con- 
ditions |P|) with respect to the measures is said to be the type-II q-Hahn 
multiple orthogonal polynomial of a multi-index n G W , i.e., the conditions 

N 

^p|' ao ' JV ( s )(4 fc] < i ' ao ' JV ( s )v^ + l) = o, k<m-i, i = i,...,r, 

(11) 
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hold, where the symbol (s)^ denotes the following polynomial of degree at most 
k in the variable x(s) 



k-l 



3=0 



(4 feI = II q \ Z i 1 = - 1) • • • x(s - fc + 1). 



Notice that (s)^ is a polynomial of degree k in x(s), and the orthogonal- 
ity conditions ^ have been written more conveniently as (fTTj) . When r = 1 
definition 12. II gives the standard orthogonality In addition, 

{1,( S )W,... S ( S )W}, 

constitutes a basis of the linear space of polynomials of degree at most k in x(s). 
Indeed, we consider this basis as the canonical one. When q goes to 1, (s)^ 
converges to (— 1) (— s)k, where (s)k denotes the usual Pochhammer symbol 
(s) fe = s(s + 1) • • • (s + k - 1), (s)o = 1. 

An important feature for P^' a °' N (s) is the existence of raising operators [3] 

Dait « ,N p «,« ,N {s) = _ g -^±^ M + ai + a ^ p S-S : , ao -l,N+l {s ^ (12) 

D a "°"" w = f — = — Vivr^ — v»;" a »' w (s) , 
\uj*- 1 ' a °- 1 -- v+1 («) 9 7 

where the multi-index eTj denotes the standard r dimensional unit vector with 

the ith entry equals 1 and otherwise, and V = f — — 7 -. Indeed, 

yx(s) 

D a uao ,N f{s) = A * 0>ai ,N( s . q) v /(s) + B a , ai ,N {s , g)/(s)j (13) 

where 

A *o^N {s . g) = g -^l [sUN + aQ _ 3 + l]q (W) 

B^ N (s; a) = L + a l]q[ N -. + !],- ^g^^ ) 

= ([iV + l] g [oi], - ar(a)[oi + ao]^ 12 ^) . (15) 



3 High order g-difference equation 

The next theorem 13 . II extends the result of [3j concerning the q-difference equa- 
tion. We will combine the lowering and raising operators to get a (r + l)-order 
g-difference equation. As a corollary, an explicit difference equation for the 
Hahn multiple orthogonal polynomials is given (see [7] ) . 

An explicit expression for the type-II monic g-Hahn multiple orthogonal 
polynomials can be obtained as a direct consequence of a summation by parts 
and formula (|12|). 
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Proposition 1. The following finite- difference analog of the Rodrigues formula 



P2 



&,a a ,N/ \ _ ( — 1) " 1 



(iy+o )ifli+nr = i "i 



'00 = 



T q (s + l)f q (N-s+l) 



IL=i(N +ao + a k + l\q)„ k 



T q {a Q +N-8 + 1) 



iir 9 ( ai + s + i) 



r g (a, +7i, + s + 1) \ 9—5— s r g (a + iV-s + l) 



r,(« + i)r,(Jv-|ri|-« + i) 



(16) 

holds, where (a\q)k = IIm=o[ a ~'~ m ]9 = ^<j( a + k)/T q (a) is the q-analogue of the 
Pochhammer symbol. 

Notice that the above Rodrigues-type formula characterizes the type-II q- 
Hahn multiple orthogonal polynomials in terms of a finite-difference property. 

Proof. Replacing (s)^ in ITT|) by the following finite-difference expression 

(fc-2) 

\[k] _ ±_ 



_VCs -I- ll[ fc+1 ] 

[k+l] q [ + )q 



the orthogonality conditions can be written in a more convenient way as follows 

N 

P r T ao > N (s)V(s + l)[ k+ %^ N ( s ) v x(s + I) = 0, 
k = 0, 1, . . .,m - 1, i = l,2, ...,r. 



s=0 



(17) 



From here, using the summation by parts one gets (|12p . Recursively using this 
raising operator gives the Rodrigues-type formula (fT6|) . □ 

In the next theorem we will need the following auxiliary lemma. 

Lemma 3.1. Let A be the following r-dimensional matrix 



( 



A = 



[111+01-02], 

[n 2 +a2-ai\ q [ n 2\ q 



1 



1 \ 

[m+ai— a r ] q 

[n 2 +a 2 — a r \ q 



\ 1 1 

\ |n r +o r -oi|, [7i r +a r -a 2 \ q 

1 



\7ii + ctj — a 



j J </ 



then the determinant of A is 

T— 1 r 



IT Y\\ ak ~ a ih\- n i - n k + ai - a k ] q 



det A = 



k=l 1=1 



[ \ ] \ "I "/ a k ] q 



k=l 1=1 



(18) 



Here for proving (|18p we will follow the operations indicated in Lemma 
2.8, p. 18]. 
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Proof. Let us proceed by column and row operations on the matrix A. Observe 
that, for k = 1, . . . , r and i = 2, . . . , r the following relation 



[OH - Qij 



(19) 



yields. 

Therefore, based on (fT!))) if denotes the fcth column of A (fc = 1, . . . , r) 
one gets 



det A = det(A 1 ,A 2 -Ai,...,A r - Ax) 

I 1 Sl,2 



IT ° fc > i 



\z=2 



• ai,r \ 

1 (22.2 • ' ' 02,r 



where a k ,i = a k ^q ^ 1 [q ' 2 1 + . Now, if one substracts the first 

row from the other ones, and takes into account that 



a-k,i— a>i,i — afc,iOi,zMfe,i 
then 



Mfc.i = [ni-n k +oti-ak] q i,k = 2,...,r, 



detA= ( [Jay J ( aij ) (H<kiMi,i 

U=2 



/ »2,2 02,3 ' ■ ' 02,r ^ 
0-3,2 ^3,3 • ■ ■ a3 jr 



\k=l 



where X h i = [an - a{\ q . 

Finally, repeating the previous column and row operations -but on lower 
dimensional matrices- the expression (|18[) can be inductively proved. □ 

Theorem 3.1. The type-II monic q-Hahn multiple orthogonal polynomial P~' a °' N (s) 
is an eigenfunction of the following (r + l)-order q-difference equation 

' A 



Ax(s) 



(20) 



i\\n\ + a + cti 



cn+1, ao+l,JV— 1 1 pS,a ,iV 



whe 



^ _ /y-» {-l) k+l [n k + a k + a + l] q ]li=l,i/i [ n * + a fc ~ 

\fc=i n[=i,i/fc [«* + a; - n fe - a fc ] 9 t nfc + ak ~ n i ~ a ^ 



1 2 IIfc=l IlLl^fc ~ a l\q\ n l - n k +Ctl - ak]c, 
Uk=lUl=l[ n l + a l - a k]q 



i = l,...,r. (21) 
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Proof. Taking into account the expressions (fT2|) - (fT5)) . the g-Hahn multiple or- 
thogonal polynomial can be expressed in term of the raising operator as follows 



pS,a ,N ^ _ Q 



\ +\H\+a -l 



L |n| + on + a ] q 

(_ B ao+l,a < +l,JV-l( fl . q ) X + A ao+l,ai+l,JV-l (s . p 3+e,, ao +l,N-l ^ 

where i = 1, 2, . . . , r. Hence, for k < r, one gets the following relation 

N 



£ P,r- JV ( s )( s )^- 1 ]^+ 1 ^^( s ) A x(.s - |) 

AT 

^(s)i" t " 11 < fc+1 ' ao,JV (s) A x(s - i) (22) 



s=0 

w + l"l+°o-i AT 



]n| + a, + a ] q s=Q 



^ a „+l,a. + l,N-l (s;?)I + A a„ + l !a , + l, W -l (s;(;)v |p5^ !ao + l, N -l (s ^ 



Now, transforming (|22|) by doing a suitable combination of the orthogonalizing 
weight v% k+1 ' ao ' N (s) with the terms A ao+1 ^' +1 ' Ar - 1 (s; g) and s^o+i^i+i^-i^. q ^ 
and summing by parts 



N 



q e [n k + a k - ai 



E^-"W(*'-V +I -"WA,( 8 -«- (|fl|+a(+ao + 1]i 

^pa+a«+i,iv-i (8)(8) K-i]„j, + i„ + i,»-i (s) A ^..j), 



s=0 

N 



X 
s=0 

(23) 

yields, where = JV +l"l+"fc+ f »fc+ a o _ jj er6j a i me d to shift the parameters ao and 
N in the orthogonalizing weight v q k+1,a °' N (s) we have used the relations 

q fc + l,ao + l, j V-l (a) 

= q2 [ N - g] qi 

Wo ' (s) 

t£*+i,"o+Mr-i( s ) ( s + [TV + a - a]. 



«g (s + 1) g 2 [ S ) q J [a fc + s + 2J 

By using recursively relation (|23p one gets 

JV 



(24) 



E ^+5 ! ' a ° +1 '^ 1 ( s )( s )i"^ 11 C +1 " 0+1 ' JV " 1 ( s ) A - 4) 

iV 

= »ME^ ao+r,,f "»( s )i"" 11 »r +I ' , ' +r ' , '"( s ) A *(» - 1)» 

where e = X/I=i ^ an d 

_ + q + a; + 1] / -j^ [q fc - aj + n fc ] Q \ _ 

° W ~ g-C-l)(«-D-i [a* - m + n k ] q I f = \ [|n| + a + a, + 1] A ' *' ~ A ' " ' ' ' " 
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Now, based on (fTO|) wc will find that there exists constants {£i}[ =1 (not all 
zero) such that the relation 

Ap | ) «o,iv (s) = |^ eiP |_ + |,ao + i,iV-i (s)) A =^, (25) 

is valid. Thus, for finding explicitly £1,. .. ,£ r one takes into account (fM|) and 
^ to get 



N 

(r \ AT 

E p l ^r o+r ' w ~>)(s)^S Qfc+1 ' QO+r ' w ~ r ( s ) Ax( S - j). 
Z=l / s=0 

(26) 

Now, the left hand side of this equation can be easily transformed by means of 
a summation by parts and orthogonality relation pip . Indeed, 

N 



[n k + a + a k + l] q 



N 



q" 2 



Based on (|23|) and (|24| the right hand side of this expression transforms into 
the equation 



N 

s=0 
N 

6 fc E p l-l' ao+r,JV_r ( s )( s )l"' t_1] < it+1 ' ao+r ' JV " r ( s ) - *)> 



(27) 



s=0 

where 



From equations ([2"B"|) and (j2"7|) leads the following linear system of equation for 
the unknown coefficients £i, . . . , £ r , 

r 

fe fc = E^ s M, k = l,...,r, S£ = 6, £ = (&,...,&•), (28) 

i=i 

where the entries of the vector b and matrix S are as follows 

bk = [rik + ctk + &o + 1] 9 , Sk,l = a>k,l[\™\ + a o + + 



By the Cramers rule, the above system (|28|) has a unique solution if and only if 
the determinant of 5 is different from zero. Observe that S = A ■ D, where D 
denotes the diagonal matrix 

D = (d k ,i)l j =1 , d k ,i = [\n\ +a + ai + l] q S k ,i. 
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Thus, from lemma I3TTI formula (fT5|) . one gets 

detS = det(A ■ D) = (f[[\n\ + a + on + 1] 9 ] det A ^ 0. 



Accordingly, if C'ij is the cofactor of the entry Sjj, and Sj(b) denotes the matrix 
obtained from S replacing its jth column by b, then 



det^) 
51 detS ' 



where 

r 

det S l {b)=Y,b k C, 



k,l 
k=l 

(-i) fc+ 'riuw^+^ 



Vi*-^ — 

fe =l n i= i,i^fc [ n i +<*i- n k - a k ] q Y\ r j=k+1 [n k + a k - nj - a 3 } q 



Consequently, expression (j2lj) yields. 

Finally, applying the following product of r operators (nl=i Z3 Qi+1 ' Qo+1 ' Ar_1 ) 
on both sides of the equation (|25|) and considering that these raising operators 
are commuting, the expression (f20|) holds. □ 



In particular when r = 2, the g-Hahn multiple orthogonal polynomial verify 
a third order g-difference equation (see [3J theorem 2.1, pp. 7-8] and next 
corollary). 

Corollary 3.2. The type- II monic q-Hahn multiple orthogonal polynomial ^"/'" 2 2:Q °'^ 
verify the following third order q-difference equation 

a 4 (s)AV 2 ?/ + a 3 (s)AVy + a 2 (s)Ay + ai(s)Vy + a (s)y = 0, (29) 

where 

a 4 (s) = q 23-3 A *o+l,*x+l,N-l( s . q - )A a +l,* 2 +l,N-l( s _ 1; g ^ 

/ Ra + l : Q2 + l:A r -l('„ _ 1 . „1 4a + l,a2 + l,A r -l(V n \ 

a 3 (s) = ^o+i,c 1+1 ^- 1(s . } ( * (j M) + ± (*,g) 

\ q L s q L s 

+ q s - 1 B ao+1 ^ +1 - N - 1 ( S ; q)A a " +1 ^ + 1 ' N - 1 (s- 1 q), 
a 2 {s) = B«o+i, Q i+i,JV-i (s . g ) B «o+i^2+i^-i( s; g ) 

w+ i fl i +ao - 2s+3 /^A^+^+^-^s;?) 6^ Q0+1 ' Ql+1,JV_1 (s;g)\ 



ai(s) = g~ 
ao(s) = 9" 



V [N + "l + + l]^ 1 [\ft\ + a 2 + a + l]^ 1 J ' 

N+N+^o+i /^ 1 B ao + 1 ' a3+1 > iV '- 1 (s;g) ^•B ao+1 ' ai+1 ' JV ~ 1 ( s ;?) > \ 
\ [|n| + ai + a + l]" 1 [|n| + a 2 + a + l]" 1 / ' 



and ^1 - q 2 ,„,_„,, , £2 - q 2 [„,_„.,!„ 
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The above g-difference equation can be considered as an extension of the 
hypergeometric-type difference equation Again here the combination of 

lowering and rising operators (|13[) was the key fact to obtain a third order 
difference equations having g-Hahn multiple orthogonal polynomials as eigen- 
functions [3]. 

Remark 3.3. Notice that from the q-difference operator i20\) can be obtained 
the difference equation studied in J^j for the Hahn multiple orthogonal polyno- 
mials since when q goes to 1 the lattice x(s) transforms into a linear one s. In 
particular, for the above third order q-difference equation \29\) one gets the same 
type of third order difference equation with the following polynomial coefficients 

a 4 (s) = s(s - l)(a +N-S+ l)(a + N - s + 2), 

a 3 (s) = s(a +N-S + 1) [2a + a 2 + 4 + (ax + a 2 + 3)N 

-(2a + ai + a 2 + 6)s] , 
a 2 (s) = [(an + l)N - (a + a x + 2)s][(a 2 + 1)N - (a Q + a 2 + 2)s] 

-(a +a 2 + 2)(a + N — s + l)s, 
ai(s) = [ni(a + a.\ + N + 1) + (a + a 2 + N + l)n 2 - nin 2 ](a + N — s + l)s, 
a (s) = N[(a 2 + l)(a + «i + N + l)n x + (a x + l)(a + a 2 + N+ l)n 2 

+ (a + N)nin 2 ] - [(a + a 2 + 2)(a + ot\ + N + l)ni 

+ (a + a x + 2)(a + a 2 + N + l)n 2 + (N - l)mn 2 ]s. 



4 Limiting transitions 

Observe that when q goes to 1 the expression (TTHJ) transforms into the Rodri- 
gues-type formula for Hahn multiple orthogonal polynomials in the linear lattice 
x(s) = s (see [1]) 

H s, ao , N( , = (-i) 1 " 1 r( s + i)r(jv - s + 1) 

" U r k =A\n\+a + a k + l) nh T(a + N-s + l) 

(ft n^hrrr^ + - + ' + ") n.lm»-to-U m 

Now, based on this Rodrigues-type formula another limiting transition be- 
tween discrete and continuous multiple orthogonal polynomials can also be ob- 
tained. The basic tool for establishing this limiting transition is the usual ap- 
proximation of derivatives by means of finite-differences. 

Suppose that f(s) is a function defined on an interval of the real line, which 
contains the linear lattice {si}fL Q . Furthermore, f(s) possesses first derivative 
on {si}fL Q , and second derivative for every x% £ [ s % ~ Sj), i = 1, . . . , N. Thus, 

vf{si) = = / [Si) - -f {Xi) 

= f'( Si ) + O(h), 

yields. In general, if f(s) has n derivatives at points and 2n derivatives for 
any \i one gets 

V"/( Si )=/ (n) (^)+0(/* n )- (31) 

Notice that the change of variable s — Nx, transforms the interval [0, N] 
into [0, 1], which is precisely the interval of orthogonality for the multiple Jacobi 
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polynomials studied in [TU1 HU [H] ■ The simultaneous orthogonality conditions 
were considered with respect to an AT-system of Jacobi weights on [0, 1] with 
different singularities at and the same singularity at 1. Below we will show 
that when TV tends to infinity, i.e., the step h = Ax = 1/N in the new variable 
x tends to 0, the Hahn multiple orthogonal polynomials (|30|) will tend to the 



aforementioned monic multiple Jacobi polynomials P~ a °' a \x) (see below the 



explicit expression ((411) ). 

Proposition 2. The following limiting relation is valid: 

lim N-WH°t' ao > N (Nx) = P*' a °(x). (32) 

N^oc 

Proof. For simplicity let us consider the multi-index n = (n>i,ri2) as well as the 
interval [0, AT — 1] as the support of the orthogonality measures (fTT)|) . The proof 
for r > 2, i.e., n — (rii, n 2 , . . . , n r ) follows the same steps described below. 

Firstly, let us show that for N large enough the term T(N — Nx)/T(ao + 
N — Nx) contained in (f5Uj) behaves like 

(1 - x)- ao N- ao . (33) 

This estimation follows immediately from the well known asymptotic relation 
for the gamma-function [9] 



r(z + a) 

r( 2 ) 



1 + 



|argz| <w-S, 5>0. (34) 



Second, for N large enough we will prove that the expression 

/ A T(Nx + l) T(a t +n l + Nx + l) \ T(N + a - Nx) 

r ( a i +Nx + 1) T(Nx + 1) J T(N - \n\ - Nx) ' 

(35) 

behaves like 

N ao (j^x~ ai ^x ai+n *\ (I-*) 004 " 1 * 11 - ( 36 ) 
Indeed, using the relation for the nth backward difference 

V"»(») - E mZT&M* - fc ) = E ^P** " (37) 
the expression (f35|) becomes 

jy-i"ir(7Vx + 1) y> y> (-n 2 )j (-m)fc r(jv(a; - fc) + ai + ni + 1) 

T(Nx + ai + 1) ^ ^ j! /fc! r(JV(a; - fc) + 1) 

r(A^(a; - fc) + l)r(JV(a; - j - fc) + a 2 + n 2 + 1)T(N + a a - N(x -j- fc)) 
T(N(x -k) + a 2 + l)r(iV(a; - j - fc) + 1)N - \n\ - N{x - j - fc)) ' 

(38) 

Hence, taking into account p4|) the above expression (|38|) can be rewritten as 

(Nx)-^ E " 3 1>k (N(x - k)) a ^ +ni (N(x - fc))"" 2 

j=ofc=o J - k - (39) 

(iv( s - j - fc)r+" 2 (jv - tv(x - j - fc)n+i"i + o (Jp V 

13 



Using again (|5T)) . but this time to express in terms of the backward differ- 
ence operators, i.e., 



N a " (jJx~ a >\7 ,H x a * +n ^ (1 - x)°"'+W + O (J^ ) . (40) 



Now, according to (|31|1 one can express the backward difference operators in- 
volved in (|40[) in terms of the ordinary derivatives. Thus, for N large enough 
one verifies that indeed ([55]) behaves like (|36p. 

Finally, from and the above estimations and the proposition 
holds. □ 

Recall that the multiple Jacobi polynomials are given explicitly by the 
Rodrigues-type formula [12] 

pa.ao/ \ _ (-l) |ra| (l - a) / TT -aj d n * ai+rli \ , _-\a +\n\ 

H [ ,_ nLi(i«i+«o +«*+!)«, IM ^ / ' 

(41) 

In fact, P~' a °(x) verifies the following orthogonality conditions 

P?' a °(x)x a >(l~x) a °x k dx = Q, k = 0,l,...,m-l, i = l,2,...,r. 

Regarding the differential equation that these polynomials satisfy we refer 
to [5]. Observe that this case constitutes a special limiting case of ([20)) . 



5 Conclusion and future directions 

To the best of our knowledge there is not in the literature any other high order 
g-difference equation different from (|20[) having multiple orthogonal polynomials 
-with g-discrete orthogonality- as eigenfunctions. Furthermore, the q-difference 
equation obtained here is quite versatile since other difference and differen- 
tial equations can be simply obtained as a limiting case. For instance, the 
difference and differential equations for the multiple Hahn and Jacobi polyno- 
mials, respectively are examples of these limiting cases. Also the well known 
hypergeometric-type difference equation for Hahn polynomials as well as the hy- 
pergeometric equation for Jacobi polynomials are particular cases when r = 1. 
However, more general situations demand our attention. Firstly, the multiple 
orthogonal polynomials on the lattice x(s) = c\q s + C2<i s + C3, where ci, C2 
and C3 are constants independent on s must be considered in the same fashion 
that here. Finally, more general systems of measures such that under certain re- 
strictions on their elements (measures) one can recover the g-difference equation 
(f2"0"f must be analyzed. In this direction, the Askey- Wilson multiple orthogonal 
polynomials could be an interesting challenge to be considered. 

In closing, this paper outlines the important points and techniques to be 
followed in future investigations aimed to deduce the high order difference equa- 
tions for g-Charlier, g-Kravchuk and g-Meixner multiple orthogonal polynomials. 
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